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Abstract. According to recent theoretical studies, the progenitors of Long Gamma Ray 
Bursts should be very fast rotating stars, massive enough but not so for collapsing into a 
black hole. In addition, recent observations seem to show that stars of about 20 solar masses 
could be at the origin of LGRBs. At low metallicity B-type stars rotate faster than at higher 
metallicity. We found with the ESO-WFI an occurrence of Be/Oe stars, that are quasi criti- 
cal rotators, 3 to 5 times larger in the SMC than in the Galaxy. According to our results, and 
using observational clues on the SMC WR stars, as well as the theoretical predictions of 
the characteristics must have the LGRB progenitors, we have identified the low metallicity 
massive Be/Oe stars as potential LGRB progenitors. To support this identification, the ex- 
pected rates and the numbers of LGRB were then calculated and compared to the observed 
ones: 3 to 6 LGRBs were found in the local universe in 11 years while 8 were actually 
observed. 
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1. Introduction 

The first section of this document deals with 
the Be phenomenon, their rates, their rotational 
velocities in low metallicity environment. The 
second section deals with the long gamma ray 



burst or type 2 bursts and their links with the 
massive Be and Oe stars at low metallicity, 
according to the last theoretical developments 
and observational facts. 
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2. Be stars in low metallicity 
environment 

The Be-phenomenon concerns the main se- 
quence OBA-type stars, which spectrum has 
displayed at least once in its life emission 
lines mainly in the hydrogen. These lines 
come from a rotationally supported circumstel- 
lar disk formed by matter ejections of the cen- 
tral star. Some information regarding the Be 
phenomenon can be found in Porter & Rivinius 
(2003). Nowadays, we know that Be stars are 
rotating very fast, very close to their critical 
rotational velocity but it is not clear whether 
an additional mechanism is needed for eject- 
ing the matter. It was also reported that the Be- 
phenomenon seems to depend on the metal- 
licity too (Maeder et al. 1999) and evolution- 
ary status (e.g., Fabregat & Torrejdn 2000; 
Martayan et al. 2007). A possibility for inves- 
tigating these effects is the study of open clus- 
ters. In the Milky Way a reference study is the 
one by McSwain & Gies (2005). 

2.1. Ratios of Be stars at low metallicity 

Maeder et al. (1999) and Wisniewski et al. 
(2007) found that the number of Be stars 
by open clusters seems to increase with 
the decrease of the environment metallicity. 
However, the number of stars and open clus- 
ters used is relatively small and it was neces- 
sary to improve the statistics by increasing the 
number of clusters, to quantify and compare 
the trends/ratios between the Small Magellanic 
Cloud (a low metallicity galaxy) and the Milky 
Way. Martayan et al. (2010) did a spectro- 
scopic survey of the SMC with the aim to find 
the emission-line stars and the Oe/Be/Ae stars. 
The ESQ WFI (Baade et al. 1999) in its slit- 
less mode was used and about 3 million spec- 
tra were obtained. The spectra of 4300 stars in 
SMC open clusters were extracted, analyzed, 
and the emission line stars found (for more de- 
tails see Martayan et al. (2008, 2010)). Then 
by cross-correlation with the OGLE catalogues 
(Udalski et al. 2008) the stars were classified in 
spectral type. With the freedom degrees con- 
strained (age, metallicity, etc) the rates of Be 
stars to B type stars per spectral-type category 



were computed in the SMC and in the Galaxy 
with data of McSwain & Gies (2005). The 
comparison of the Be to B stars per spectral 
type category rates is shown in Fig. 1. 
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Fig. 1. Ratios of Be stars to all B-type stars per 
spectral type categories. The SMC sample is only 
complete till B3 spectral type. The SMC ratios are 
compared to those from McSwain & Gies (2005) in 
the Milky Way. 



The SMC sample is complete till spectral 
type B3. For spectral type ranging from BO to 
B3 the rates in Fig. 1 indicate that there are 3 
to 5 times more Be stars in low metallicity en- 
vironment than in our Galaxy. 

2.2. rotational velocities at low metallicity 

As shown in Sect. 2. 1, the occurrence of the Be 
phenomenon is larger in lower metallicity envi- 
ronment. It is explained by the fact that at lower 
metallicity O, B, A, and Be stars rotate faster 
than their counterparts in higher metallicity en- 
vironment (Keller 2004; Martayan et al. 2006, 
2007; Hunter et al. 2008). It is due to their 
lower mass-loss at low metallicity (Bouret et 
al. 2003; Vink 2007) that implies a lower angu- 
lar momentum loss (Maeder & Meynet 2001). 

With FLAMES observations Martayan et 
al. (2007) determined the rotational velocities 
of SMC Be, B stars. They also determined their 
ZAMS rotational velocities taking into account 
the fast rotation effects such as the gravita- 
tional darkening (Fremat et al. 2005; Zorec 
et al. 2005). The ZAMS rotational velocities 
of SMC Be stars are compared with theoreti- 
cal models from Ekstrom et al. (2008) in Fig. 



Martayan et al.: Massive SMC Be/Oe stars and LGRBs 



3 



2. Their models seem to be able to reproduce 
properly the observations. 

One can also note that SMC Be stars at 
their birth rotate very fast, with ZAMS rota- 
tional velocities exceeding 550 kms 'for the 
most massive of them. However, for such very 
fast rotators, the stellar evolution could be 
modified and instead of following a classi- 
cal stellar evolution they could follow a quasi 
chemically homogeneous evolutionary as de- 
scribed in Brott et al. (2011). In Fig. 2, a box 
with dashed lines shows the predicted area of 
stars following this quasi chemically homoge- 
neous evolution. 
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Fig. 2. ZAMS rotational velocities of SMC Be and 
Oe stars compared to theoretical curves from the 
models of Ekstrom et al. (2008). In addition a box 
shows the location of stars that according to the the- 
ory (Brott et al. 2011) can follow the quasi chemi- 
cally homogeneous evolution. 

Therefore one can conclude that the mas- 
sive Be and early Oe in low metallicity en- 
vironment such as the SMC could follow the 
chemically homogeneous evolution. It is worth 
to notice that Martins et al. (2009) found SMC 
WR with properties that could be explained 
only from the chemically homogeneous evolu- 
tion. 

3. LGRBs and Oe/Be stars 

This section deals with the type 2 burst also 
called long soft gamma ray burst (LGRB) 
because they are longer than l-3s. Woosley 
(1993) proposed to explain the LGRB, the 
model of a massive fast rotating star collaps- 
ing into a black-hole during SNIb,c explosion. 



The energy for the explosion is provided by the 
potential energy of the matter falling onto the 
collapsed stellar core. This remaining energy 
will not accrete onto the black hole that implies 
for the infalling matter to have enough angular 
momentum to remain in a disk before accre- 
tion. 

3.1. Observational and theoretical facts 

Iwamoto et al. (1998, 2000) support the idea 
that massive fast rotating stars are at the ori- 
gin of the LGRBs. Thone et al. (2008) found 
that the LGRB GRB060505 is hosted in a low- 
metallicity galaxy, with a high star-formation 
rate. The LGRB came from a young environ- 
ment (6 Myears) and from an object having 
about 32 Mq. Campana et al. (2008) found 
that the LGRB0602 1 8 progenitor had an initial 
mass of 20 Mq. 

Theoretically, the rotation seems to be a 
key point to understand the appearance of 
GRBs (Woosley 1993; Hirschi et al. 2005; 
Fryer et al. 2007). To keep a large amount of 
angular momentum up to the last evolution- 
ary phases before the collapse, GRBs progeni- 
tors should be massive objects with low initial 
metallicities. According to Yoon et al. (2006) 
WR stars with metallicities Z < 0.002 can be 
progenitors of GRBs. 

Due to fast rotation and the rotational mix- 
ing of chemical elements, massive stars can 
undergo quasi-chemically homogeneous evo- 
lution to end up as helium WR stars satisfy- 
ing the requirements for the collapsar scenario 
(Yoon et al. 2006). Yoon et al. (2006) have 
calculated diagrams of LGRBs progenitors at 
different metallicities (including the SMC one) 
as a function of their ZAMS rotational veloc- 
ities and masses for magnetized massive stars 
following the quasi-chemically homogeneous 
evolution. 

From their models, at low metallicity the 
WR phenomenon can appear in stars having 
lower masses than those in the Milky Way. We 
recall the observations by Martins et al. (2009) 
of several SMC WR stars, whose evolution- 
ary status and chemical properties can be un- 
derstood if they are fast rotators following this 
quasi-chemically homogeneous evolution. 
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According to the different criteria shown 
by the observations and the theory, the stars 
must be massive or of intermediate mass, fast 
rotator, formed in low metallicity environ- 
ments. Complying with these criteria and the 
diagrams of LGRB progenitors massive SMC 
Be and Oe stars could become WR star and 
give a LGRB. 

3.2. Predicted rates of LGRBs from 
Oe/Be populations 

For testing the hypothesis of the SMC massive 
Be and Oe stars as LGRB progenitors, one can 
predict the rates of LGRBs based on that spe- 
cific stellar population. One can use the SMC 
as a test-bed galaxy, the number of OB stars is 
determined with the OGLE catalogues. Then 
the number of the SMC Be/Oe stars is deter- 
mined through the rates shown in Sect. 2. 1, for 
more details, see also Martayan et al. (2010). 

A base rate of LGRB/year is then obtained. 
However the gamma rays are strongly coUi- 
mated. One has to take into account the beam- 
ing angle to compute the probabiUty to see a 
GRB. 

The table 1 gives the LGRB base rates and 
the LGRB rates depending on the beaming an- 
gle considered and for different mass category 
of stars (from B2e to 08e). According to the 
mass calibration of Huang & Gies (2006) the 
stars of spectral type BO and higher (masses 
above 20 Mq) should be considered here. 

However, one can take into account the 
LGRB phenomenon beaming angle distribu- 
tion (Watson et al. 2006) instead of discrete an- 
gle values. In such case, one can get: A^l^rb ~ 
(2 - 5) X 10"^ LGRBs/yr/galaxy that can 
be compared to the observed rate of LGRB: 
^LGRB ~ (0.2 - 3) X 10-^ LGRBs/yr/galaxy 
(Podsiadlowski et al. 2004; Zhang & Meszaros 
2004). 

3.3. Predicted number of LGRBs in tfie 
iocai universe 

Another possibility is to predict the num- 
ber of LGRBs in the local universe, i. e, at 
redshift<0.2 for which the galaxy catalogue 



of Skrutskie et al. (2006) is complete. Up to 
redshift 0.5 the number of Im is about 17% 
(Rocca-Volmerange et al. 2007). 

One can get the predicted number of 
LGRBs in the local universe during 1 1 years 
(between 1998 and 2008) using the previous 
predicted rates of LGRBs. The result is shown 
in table 2 for numbers by mass category and 
discrete values of beaming angles. 

Ditto, one can also use the beaming angle 
distribution, in such case the predicted num- 
bers of LGRBs are: N^^^ -3-6 LGRBs 
(for mass categories starting from BOe to 08e). 
This number could be compared to the ob- 
served number of LGRBs in the local uni- 
verse in 11 years with data of the the GRBox 
from the University of Cahfornia at Berkeley' 
for years between 1998 and 2008, there are 8 
LGRBs. 

Therefore there is a relatively good agree- 
ment between our predictions and the observed 
rates/numbers of LGRBs in the local universe. 

3.4. Binaries 

Malesani et al. (this volume) report that sev- 
eral LGRBs occasionally occurred in "high" 
metalhcity environment. If confirmed it indi- 
cates that the usual models of LGRBs, which 
need low metallicity environment cannot han- 
dle them. Therefore, the binary models of 
LGRB creation (Cantiello et al. 2007) could 
be of interest too. Considering too the binary 
system as potential progenitor of LGRB, one 
could multiply the rates we estimated above by 
1/0.7, which gives a number of 9 LGRBs in 
the local universe (to be compared with the 8 
LGRBs actually observed). 

3.5. Other biases 

Other possible biases than the binaries could be 
taken into account such as the missed GRBs by 
the satelhtes, the obscured GRBs not detected 
as GRB, etc. On the other side, one could also 
take into account the effect of the very fast 
rotating non emission line star, the so called 
Bn stars on the statistics. For more details, the 

' see http://lyra.berkeley.edu/grbox/grbox.php 
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reader could consult the paper by Martayan et 
al. (2010). 
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Table 1. Predicted LGRB rates 



Mass category 


LGRB base rate 


Proba 5deg 




Proba lOdeg 


Proba 15deg 


B2e to 08e 


3.9-4.8 xlO-'* 


2.5-3.1 xlO" 


-7 


1.0-1.2 xlO-' 


2.3-2.8 xlO"* 


BletoOSe 


2.6-3.0x10'' 


1.7-1.9 xlO 


-7 


6.7-7.7x10-^ 


1.5-1.7 xlO'** 


BOe to 08e 


1.7-1.8 xlO"" 


1.1-1.2x10" 


-7 


4.4-4.6 xlO-'' 


0.98-1.0 xlO-« 


09e to 08e 


5.3-6.2x10-5 


3.4-4.0 xlO- 


-8 


1.4-1.6x10'^ 


3.1-3.6x10-' 


OSe 


2.4-2.8 xlO-5 


1.6-1.8 xlO- 


-8 


6.2-7.2 xlO-* 


1.4-1.6 xlO-' 



Table 2. Predicted LGRB numbers in 1 1 years in the local universe (redshift<0.2). 



Mass category 


Niunber for angle=5'' 


Number for angle=10° 


Number for angle=15° 


B2e to 08e 


3-4 


11-14 


25-31 


Bleto08e 


2-2 


7-9 


16-19 


BOe to 08e 


1-1 


5-5 


11-11 


09e to 08e 


0-0 


2-2 


3-4 


08e 


0-0 


1-1 


2-2 



